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I . INTRODUCTION 

Double-hctero-junction  GaAi’As/GaAs  lasers  constitute  a key 
component  in  integrated  and  fiber  optical  systems.  The  basic  configu- 
ration is  illustrated  in  Kig.  (la).  When  a strong  dc  bias  is  applied 
between  the  stripe  electrodes,  the  segment  of  the  GaAs  layer 
immediately  below  the  stripe  becomes  active,  resulting  in  lasing 
action.  The  GaAs  layer  has  a higher  dielectric  constant  than  the 
GaAMs  layer.  Hence,  it  is  easy  to  explain  the  existence  of  a guiding 
mechanism  in  the  vertical  direction  (i.e.  perpendicular  to  the  layers), 
by  means  of  a conventional  slab  waveguide  analysis  [2],  However, 
it  has  been  experimentally  observed  that  the  field  is  also  confined  in 
the  transverse  direction  (i.e.,  parallel  to  the  layers).  This  effect 
cannot  be  explained  in  terms  of  an  infinite  slab  waveguide  since 
no  physical  mechanism  is  present  for  transverse  field  confinement. 

it  is  now  generally  accepted  that  the  dielectric  properties 

of  the  active  region  are  slightly  altered  by  the  lasing  action.  It 

has  been  assumed  by  some  workers  that  lasing  action  causes  a slight 

increase  of  the  relative  dielectric  constant  in  the  active  region 
-3 

(of  the  order  of  10  ),  which  is  sufficient  to  explain  field  confinement. 

Only  the  real  part  of  the  dielectric  constant  was  considered  in  explaining 
transverse  field  confinement.  However,  some  recent  investigations  [3], 
indicate  the  real  part  of  the  dielectric  constant  in  the  active  region 
need  not  be  higher  and,  indeed,  may  even  be  lower  than  that  in  the 
passive  region.  Modal  confinement  is  shown  to  be  possible  by  the  gain 


action  itself. 
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In  previous  analysis  [3],  the  x and  y directions  were  treated 
independently.  For  instance,  in  [3],  only  the  structure  which  is  infinite 
in  extent  in  the  y direction  is  considered.  In  the  analysis  in  this 
paper,  however,  the  finiteness  in  the  y direction  is  taken  into  account  by 
the  use  of  the  concept  of  effective  dielectric  constants  (EDC),  which  has 
been  successfully  used  to  analyze  a number  of  lossless  passive  dielectric 
waveguides  and  components  [4] , [5] . 

In  the  following  sections  we  will  describe  the  analysis 
of  the  laser  structure  and  a number  of  numerical  data  for  the 


fundamental  mode  will  be  presented. 


I 
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II.  METHOD  OF  ANALYSIS 

The  basic  geometry  of  the  GaAs  laser  is  shown  in  Fig.  (la). 

A GaAs  layer  of  thickness  2d  is  sandwiched  between  two  identical 
GaAfAs  layers  of  thickness  h.  When  an  appropriate  d.c  bias  is 
applied  between  the  electrodes  A and  A' , the  GaAs  region  under  the 
electrodes  becomes  active  and  causes  a laser  action.  We  model  this 
effect  by  assigning  the  region  under  the  electrode  a complex  permit- 
tivity ej  + je^,  as  shown  in  Fig.  (lb).  In  this  figure  the  presence 
of  the  electrodes  is  ignored  as  its  effect  on  the  optical  field  is 
negligible.  From  Fig.  (la)  we  observe  that  the  field  is  symmetric  with 
respect  to  the  y axis  and  antisymmetric  with  respect  to  the  x axis. 
rherefore  a magnetic  and  electric  wall  can  be  placed  in  the  y and  x 
planes,  respectively,  without  changing  the  field  configurations. 

It  is,  hence,  only  necessary  to  consider  the  quadrant  illustrated  in 
Fig.  (2).  Along  the  x direction  we  distinguish  two  regions:  I.  which 

contains  the  active  material  and  II,  which  is  purely  passive. 

For  the  present,  all  dielectrics  are  considered  lossless. 

Region  I contains  three  dielectrics:  (1)  the  active  GaAs 

with  complex  dielectric  constant 

ei  = ei  + J el  (la) 

0 < £'i  <<:  el  (lb) 

where  e > 0 represents  the  gain  mechamism. 

(2)  the  lossless  GaAiAs  with  real  > 0 and 

(3)  the  air  with  = 1.  Also  we  assume  the  following  relations  to  hold: 

I 

I 


\V 


Fig.  lb:  Waveguide  model  of  the  GaAs  laser. 
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MAGNETIC 
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e'2  “ c2 

(2a) 

e;>e2 

(2b) 

el  > *2 

(2c) 

A dielectric  waveguide,  as  shown  in  Fig.  (lb)  is  known  to 


support  the  propagation  of  hybrid  modes,  classified  into  two  possible 


x Y 

field  configurations:  E and  E . 

6 H 

two  scalar  potentials,  0 and  0 . 


These  modes  can  be  represented  by 
-jk  z 

Assuming  a e z dependence, 


Maxwells  equations  [1]  in  terms  of  these  potentials  are: 

E . i^!+wk0» 
x oyox  z 


(3a) 


H = (k 

y z 


si,*" 

3x2 


(3b) 


c = 5#.  ,(nu  $£ 

Ez  e 3y  J dx 

r J 


(3c) 


E = f <kz  ' H)0e 

y e Z - l 

} r 3x 


(4a) 


H = - 'ue  k 0e  + 
x o z oyax 


(4b) 


where 


H 

z 


30e  ..  A0H 

jaJeo  • Jkz  W 


(4c) 


e = permittivity  of  free  space 
o 

= relative  permittivity  in  the  region  of 
application 


j = 


8 


However,  due  to  the  excitation  the  dominant  electric  field  is  in  the 
y direction.  We  therefore  assume  0**  = 0.  Furthermore,  because  the 
differences  between  €”1  and  as  well  as  between  and  e.^  are  small, 
we  assume  the  transverse  wave  numbers  are  small.  This  implies  that 
second  derivatives  can  be  neglected  and  the  dominant  mode  is  of  the 

y 

E type  with  the  principal  field  components  E , H^,  and  E 

A rigorous  solution  of  Maxwell's  equations  for  the  present 
laser  geometry  would  be  exceedingly  complex.  However,  it  is  possible 
to  introduce  a simplification  by  the  use  of  the  concept  of  effective 
dielectric  constant.  If  regions  I and  II  were  infinitely  wide,  they 
would  reduce  to  the  double  layered  slab  structures  as  shown  in 
Fig.  (3a)  and  Fig.  (3b),  respectively.  The  propagation  constants  for 
these  double  slab  waveguides  can  be  determined  by  matching  the 
tangential  fields  across  each  boundary.  Both  these  structures  can  then 
be  replaced  by  infinite,  homogeneous  regions  having  a effective 
dielectric  constant,  which  may  be  thought  of  as  the  dielectric  constant 
of  the  hypothetical  medium  in  which  the  phase  velocity  is  identical 
to  that  of  the  surface  wave  in  the  original  (slab)  structure.  Referring 
to  Fig.  (2),  we  replace  region  I and  II  with  vertical  slabs  of  effective 
dielectric  constant  ^ and  > even  though  these  regions  are  not 

infinite  in  the  x direction.  The  result  is  the  structure  shown  in 
Fig.  (3c).  We  can  now  solve  the  eigenvalue  equation  of  this  structure 
for  the  propagation  constant,  which  is  assumed  to  be  that  of  the  original 
structure.  Notice  that  since  is  complex,  must  also  be  complex 
to  represent  the  gain  mechanism. 


r 
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III.  DERIVATION  OF  THE  EIGENVALUE  EQUATIONS 

In  this  section  we  derive  the  eigenvalue  equations  for  the 

structures  shown  in  Fig.  (3a),  (3b),  and  (3c)  and  define  the  effective 

dielectric  constant  for  regions  I and  II. 

Y H 

Noting  that  for  the  E mode,  0 = 0,  we  see  from  equations 

(3a)  through  (4c)  that  the  tangential  components  to  be  matched  are 
and  Ez • The  relationships  between  and  Ez  and  0 are 


(5a) 


E ± ^ 
z ^ e dy 
r J 


(5b) 


Considering  the  structure  in  Fig.  (3b)  and  noting  that  since  E^  = 0 at 

0 

y = 0,  we  can  choose  the  following  function  for  0 . 


A cos  k „y 

y2 


d > y > 0 


<t)'~  = BCcos  h (Tl^  (y-d))  + BSsin  h (T]^  (y-d))  d+h  < y < d (6) 


JI  exp[-?2(y-d-h)] 


y ;>  d + h 

Since  the  fields  must  match  for  all  z we  also  have, 


, 2 „ i 2 ,2  . 2 ^2  _ , 2 , =-2 

k = e k - k = e.  k + T\  = k + ?>„ 
z2  lo  y2  2o2  o 2 


(7) 


where  r2  ’■s  reai  and  greater  than  zero,  k^  is  real  and  is  purely 
real  or  purely  imaginary.  Using  equation  (6)  we  match  and  E^  at 
y = d and  obtain 


A cos  k „ d = B 

y2 


- 7—  k sin  k _ d = — BS 

€l  y2  y2  k2 


(8a) 

(8b) 


1 


r 
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Similarly  at  y = d + h we  obtain 


B^oshC^h)  + BSsinh(Tl2h)  = C 


— B sinh(TLh)  + — B cosh(TLh)  = - ? C (9b) 

G 2 ^ ^ ^ 

After  some  manipulation  we  obtain  the  eigenvalue  equation 

e TL  tan h (TLh)  + e??2 

k , tan  k _ d = T|  — =■  — — -f  ■-  (10) 

y2  y2  2 e2  c2  f tan  h (T^h)  +T12 

Together  with  equation  (7)  we  can  solve  for  k^2  and  define  the  effective 
dielecLric  constant  for  Fig.  (3b)  as 


/ y2 

" ei  ‘ 


A similar  analysis  for  Figure  (3a)  leads  to  the  eigenvalue  equation 

e T)  tanhCH  h)  + e 

k . tank  , d = TL  ,r;  . n~.  "ti'~  0 

yl  yl  1 «2  e2  tanh  (Mjh)  +11^ 

Again 

, 2 - , 2 ,2  ~ ,, 2 v2  _ 2 -2 
k . s s,  k • k , * «,  k + l - k +5,  (■ 

zl  lo  yl  2o  1 o 1 

and  we  define  the  effective  dielectric  constant  for  this  structure  as 


sel  " *1  ‘ Vk  > 
o 


kyl  . 2 


Note,  however,  that  in  equations  (12)  through  (14)  k^,  kyl> 
T]j,  Fj  and  are  complex  quantities. 

However,  as  a consequence  of  (lb),  the  imaginary  parts  of  all 
quantities  are  much  smaller  than  the  corresponding  real  parts.  Writing 
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k - = k + jk 

z 1 z J z 


k i = k i + jk  i 

yl  yl  J yl 


T)  = T)'  + ill" 

1 1 J 1 


e = - + ir 

n i J i 


ei  = + j<i 


(15a) 
(15b) 
(15c) 
( 15d) 
( 15e ) 


we  take  a Taylors  series  expansion  of  both  sides  of  equation  (12)  about 
(k  j > Ignoring  second  order  and  higher  order  terms  and  equating 


real  and  imaginary  parts,  we  get 


k n tan(k  , d ) = T| 
yl  yl  1 


_ , e_[  71  j tanh(T]|h)  SL 


®2  e2  ^ tan  b (^j  h ) +Tlj 


(16) 


and 


ky  1 ! Sin  (ky  1 d ) + kyl  d c°s(ky1d)]  = cos  (k'yld) 

■ c[  sin<k^d>  ' kyl  d]  ' r7  + 71 1 ei  cos  ^ky i d>  • 77 


+ n;e;Cos(k;1d).  (17) 


where 


n'  = T]|  tan  h (T)|  h)  + 
d'  = F1  ^2  tanh(Tljh)  + e2  T)| 


(18a) 

(18b) 
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N"  = ^[tanh(lljh)  + T]|  h sec  h2  (TlJ  h)  ] + e2  ?" 

D"  = * J e2  n'1'  h sec  h2(1lj  h ) + e2  tan  h (T1^  h)  + T^e  2 
The  real  and  imaginary  parts  of  (13)  are: 


(18c) 

(18d) 


e!  k2  - k'2  = r k2  + 1112 
1 o yl  2 o 1 

- k2  + §'2 

O 1 

(19) 

c','  k2  - 2k'  - k , = 2T|  T)" 
1 o yl  yl  11 

= 9 P 1 f" 

Z5in 

(20) 

Observe  that  (15)  and  (19)  constitute  three  equations  for  the  three 
unknowns  k p 'i»  Once  these  unknowns  have  been  solved,  their 

value  can  be  substituted  into  (20)  and  (17).  Note  that  equation  (17) 
is  a linear  equation  in  k^,  as  indeed,  it  should,  since  only  the  first 
term  of  the  Taylor's  series  was  kept.  Ignoring  higher  order  terms, 
equation  (14)  may  be  written  as 


el  = cel  + jeel  = *1 


fyV* 

{—)  + ^e 

o 


1 


2kyik-yi 


(21) 


Having  defined  the  E.D.C.  in  equations  (14)  and  (21)  for  the  equivalent 
slab  structure  (Fig.  3c)  we  can  now  proceed  solving  for  the  wavenumbers  in 
the  x direction  with  a analysis  similar  to  that  in  [3].  The  tangential 
fields  we  are  interested  in  matching  now  are  and  H_^.  From  equations 
(4a)  and  (4c)  we  note 


e ~ <tr 

y 


(22) 


(23) 
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Writing  0 as 


(a  cos(u  — ) 0 < x «r  a 


0 = 


I- 


W 


B exp[-  — (x  - a)]  x >a 


(24) 


we  match  the  fields  at  x = a and  obtain 


A cos(u)  = B 


- — A sin(u)  = - — B 
a a 


(25) 

(26) 


from  which  one  gets  the  eigenvalue  equation 


w = u tan(u) 


(27) 


As  before,  by  matching  the  fields  along  the  z direction  we  also  get 
another  equation  for  u and  w 


k2  = e 
k2  €el 


, • v 2 

k2  - in; 

o \a  / 


,2  w 
e 0 k + 
e2  o Va  / 


(28) 


Notice  that  u = u^_  + ju^ , w = wf  + jw^  are  in  general  complex  numbers. 

In  the  following  sections  three  cases  will  be  considered: 

(i)  Re€gl  = ee2,  (ii)  Reeel  > eg2 , and  (iii)  Reeel  < eg2 . 
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IV.  RESULTS  AND  DISCUSSION 

In  Fig.  (4)  and  Fig.  (5)  the  solution  for  the  equivalent 
dielectric  constants  defined  in  equations  (14)  and  (21)  are  displayed  for 
different  values  of  ej  and  e^.  Since  equations  (12)  and  (16)  are 
identical,  Fig.  (4)  may  be  used  to  compute  both  e'^  and  r ^ . The  equations 
were  solved  by  iteration  using  Muller's  Method,  convergence  was  quite 
rapid  and  fairly  insensitive  to  the  starting  point.  The  program  is 
included  in  the  appendix. 

We  now  proceed  to  examine  the  first  case  (i.e.,  Res  , = e ). 

el  e2 

First  we  note  that  from  equation  (28)  we  can  solve  for 

w = V k.Ki  " « 9J  - u2}  (29) 

o el  ez 

Since  Re  equation  (24)  can  be  rewritten  as 

w = { a ^ (ImAe  ) -u^]  (30) 

o e 

where  As  Defining  v = a k y ImAs  , we  see  from  equation  (27) 

eele2  o e 

that  all  physical  dependence  of  the  transverse  wave  numbers  can  be  defined 
in  terms  of  a single  variable  |v|,  called  the  normalized  gain.  Equation 
(27)  is  then  solved,  again  using  Muller's  Method. 

Figure  (6)  and  (7)  show  the  solutions  for  w^_  and  w^  vs.  |v|. 

Figure  (6),  which  indicates  the  rate  of  field  decay  in  Region  II,  shows  that 
gain  (i.e.,  ImAe^  > 0)  induced  modes  exist,  and  are  in  fact  guided  throughout. 
That  is,  the  gain  induced  by  the  imaginary  part  of  the  effective  dielectric 
constant  is  solely  responsible  for  the  field  confinement.  From  figure  (8) 
as  well  as  from  figure  (9),  which  shows  the  percent  of  the  cross-sectional 


(or  e 


I 
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power  in  Region  I,  it  is  clear  that  the  greater  the  gain,  the  more  field 

is  concentrated  in  Region  I (active  region).  The  values  of  the  field 

intensity  at  x = a (Fig.  8)  have  been  compared  to  the  results  obtained  by 

Schlosser  [3]  and  the  agreement  is  good. 

In  Case  (ii)  the  gain  as  well  as  the  increase  in  the  real  part 

of  the  FDC  in  Region  I contributes  the  confinement  mechanism.  However, 

this  case  will  not  be  discussed  further,  because  the  confinement  due 

to  the  real  part  is  understood  as  the  surface  wave  mechanism. 

Case  (iii ) requires  more  careful  examination.  It  is  known  that 

a quasi-mode  (leaky  wave  mode)  with  low  loss  can  exist  in  a passive 

channel  waveguide,  where  the  core  material  has  lower  refractive  index  than 

the  sui rounding  medium,  provided  the  core  size  is  appreciably  larger 

than  the  wavelength  [6] . The  latter  is  always  the  case  in  the  configura- 

2a 

tions  under  study,  as  we  have  typically  — ~ 20.  As  the  epithet  leaky 
suggest,  propagation  in  such  a "mode"  is  intrinsically  accompanied  by 
power  loss  due  to  leakage  from  the  core  into  the  surrounding  medium. 

Hence,  unless  energy  is  continuously  supplied  as  the  wave  propagates,  the 
leaky  mode  will  vanish  below  a detectable  level  after  propagating  over  a 
finite  distance. 

In  the  present  situation  the  supply  mechanism  is  provided  by  gain 

in  the  region  where  wave  is  propagating.  Hence,  sufficient  gain  will 

sustain  a stable  cutoff  mode,  i.e.,  mode  exhibiting  no  gain  or  attenuation 

outside  the  core.  When  this  occurs,  w^  = 0,  as  the  field  amplitude  in 

Region  II  is  expressed  as  expf-w  (~  - 1)].  We  recall  that  the  attenuation 

r 3 

constant  of  a leaky  mode  is  approximately  inversely  proportional  to  the 
index  depression  (see  for  instance,  1.6-34  of  [6]).  Hence,  the  cutoff 

| 
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point  is  achieved  with  higher  gain  for  decreasing  refractive  index 
depression,  involving  larger  leakage. 

For  a given  refractive  index  depression,  as  the  gain  increases 
above  the  value  corresponding  to  cutoff  (i.e.,  the  amount  of  gain  needed 
to  compensate  for  leakage  losses),  then  increasing  field  confinement 
takes  place. 

On  the  other  hand,  physical  intuition  suggests  that  field  con- 
finement is  more  sensitive  to  gain  increases  for  smaller  depressions. 
This  is  easily  shown  to  be  the  case  by  means  of  the  approximate  analysis 
similar  to  that  of  [3].  We  begin  by  again  defining  v = a kQ  J ^ 


v = a k Re  Ac  1 1 - j 
o e \ 


ImAe 

— V 

ReA? 

e ' 


Since  ReA?  < 0 


ImAe  \ 1 

v-Jak„VM..I  0-J  ^ 

1mA? 

I f — <<  1 , then 

UeAe  | 
e 


Let 


v 


ImAe 
e__ 

UeAe  | 
e 


v . = a k ,/  | ReAe  I 
l o e 

ImAe 

v = a k — r 

r 2 | ReAee | 

V = v + jv. 
r l 


(31) 


(32) 


(33) 


(34) 

(35) 


then 


(36) 
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In  the  case  under  study 


there  fore , 


v 


and 


w 


r 


\ 


ImAe 
e_ 

|ReA£  | 
e 


(37) 


and  — « (38) 

d(ImA-:c)  I ReAc  | 

In  other  '.wriis  , the  slope  of  decreases  for  increasing  depressions. 

The  occurrence  of  the  above  two  effects  on  the  numerical  solution  is 

shown  in  figure  (10).  In  particular,  the  approximate  linearty  of  with 

ImAc^  implies  the  existence  of  a cross  over  point  for  two  different  values 

of  Ro£g  . Figure  (11)  shows  the  solution  for  the  imaginary  part  of  w. 

In  figure  (12)  we  show  the  field  intensity  at  x = a,  with  ImAe  as 

e 

the  variable  and  ReAe^  as  the  parameter.  From  figure  (13),  the  numerical 
results  indicate  that  for  depressions  less  than  .001  more  gain  is  required 
to  maintain  guiding  than  for  smaller  depressions.  However,  for  depressions 
larger  than  approximately  .001  the  pattern  reverses:  the  greater  the 

depression  the  less  gain  is  required  to  insure  guidance.  This  phenomena, 
just  discussed  above,  can  also  be  verified  by  a approximate  solution  of 
equa tion  (27). 

We  start  by  assuming  uf  is  close  to  tt/2  . 


where 


u = (j  f-a)  + jb 

| a | and  | b | « 1 . 


I -I  2 « 

v 1 


(39) 


Since 
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2 - 1 2 

w = v(  l - (i)  j & v l - ^ (!±)  ) 

V 'v  / \ 2 v / 


and  (27)  becomes 


/ r + a + jb 


(Vr  + lVi  )'  1 ‘ 2 l ~v — +w — J ' = (2  + a + ib)tan(j  + a + jb) 
r i 


Since  « 1 and  >1,  v^_  may  be  ignored.  Furthermore,  we  note 
that  tan(— +a+jb)  ^ -1/ (a  + jb) , and  if  |ReAe  [ is  large  enough. 


1 2+a+ib 

such  that  — , : , 


2 \ j v. 


J I «1,  equation  (39)  becomes: 


(j  + a + jb) 


1 V . = - — 

i a + jb 


Solving  now  for  a and  b we  get 


n . i 

3 ^ 1 + V . 2 


_ TT  ( Vi  \ 

b - 2 '77 

1 + V. 

1 


To  show  the  behavior  at  cutoff,  we  expand  equation  (28). 


w“  - wt  + 2jwrWj  = (a  koTReAee  + j (a  kQ)^  ImAe 


2 x 2 o, 

- u + u.  ■ 2ju  u, 
r i r i 

At  cutoff,  by  definition, w ^ = 0.  It  follows,  then,  that 


(a  k ) ImAe  = 2u  u. 
o e r i 


Using  the  results  from  equations  (41)  and  (42) 


ReA€e  =-.01 


ReAee  =“006 
ReAe,  =-,004 


Im  Aee  x 10 


Fig.  11:  Imaginary  part  of  W vs . ImAc  as  parameter. 


ImAee  [2a/X  =20] 


HP  - 10? 


Fig.  12:  Intensity  of  the  field  at  x = a vs . ImAee  with  ReAe 


as  parameter. 
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3 

imAGe  (at  cutoff)  = 2 (rr ) [ ~'\A  <47> 

o (1  +v  . ) 
i 

1 t can  be  immediately  seen  that,  since  v . > 1, 

l 

ImAe  (at  cutoff)  “ — ^ 1/./ 1 ReAe  ( . That  is,  the  larger  the  depression, 
e vi  e 

the  smaller  the  gain  necessary  to  insure  a guided  mode.  Equation  (45) 
is  plotted  along  side  the  numerical  results  in  figure  (13). 

Looking  at  equation  (32),  it  is  possible  to  define  a more 
general  set  of  variables. 

ImAe 

^7  (48) 

ReAe 

e 

and  v = — J |ReAeg|  . (49) 

Data,  similar  to  that  in  figure  (10)  through  (13),  is  presented  in  terms  of 
these  variables  in  figures  (14)  through  (17). 


T 


T 


T 


T 


Fig. 


i I i I i 

0.006  0.008  0.0] 

Re  | Ac e | [2a/X=20]  hp-i<x 


13:  ImAe  at  cutoff  vs.  ReUe 
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V.  CONCLUSIONS 

Two  dimensional  analysis,  based  on  the  concept  of  effective 
dielectric  constant,  was  presented.  Several  numerical  results  are 
provided  for  the  fundamental  mode,  illustrating  gain  induced  confinement 
in  the  transverse  direction.  Gain  induced  modal  confinement  is  shown 
possible  even  for  a leaky  wave  structure  and  cutoff  behavior  is 
investigated . 
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In  this  section  a program  to  solve  equations  (16)  and  (17)  is 
presented.  The  algorithm  used  was  Miller's  method.  Convergence 
was  quite  rapid  and  fairly  insensative  to  the  starting  point. 


this  program  solves  equations  n*>)  and  cm  using  muller's 

me  T hO.)  . *L  l u • i a m t I T I h S AWE  NORMALIZED  WITH  RESPECT  TO 
<0.  HAVELS'  GTM(WAVFL)  ,M,0,E?  MUST  RE  SUPPLIED,  I*AX  IS  THE 
max  I MUM  n,jm„e<,  op  ITERATIONS.  A INITIAL  GjESS  Z3  MUST  BE 
SUPPLIED  *S  well  AS  T mE  EXPECTEO  ERROR  (ERROR)  OF  THE 
G’JE  SS . 

SI  INOICATFS  WHETHER  ETa'  is  RFAL  OF  IMAGINARY, 

E p 3 1 AND  EPS?  AWE  CRITERIA  FOR  STOPPING  ThE  ITERATION 


COMMON  E2,El ,H,0,KAVEL,PI,Sl 

pI*i.l«lS 

E?«l?.32 

PII*?.»PI 

WAVEL  ej .E-6 

h*S. »wAVFl 

D« ,?S*ha V£L 

C ncpmalIZE  with  RESPECT  TO  KB 

H*fFTI/ftAvFl.  )«H 
U*fPII/-AvEL)«0 
EPS1«1 ,E-5 
E°S2»i ,E-9 
IM,aa»S0 
ERROR*. 911 

c t no  of  data  block 


c BEGIN  root  FINDING  ROUTINE 

C TmF  fqijaTIOn  hill  be  SOLVED  FOR  DIFFERENT  VALUES  OF  El 

n9  IDO  K*1 , 1 
El*l?,T-C<-l)*,Bl 
NC • UC ♦ 1 
NC2*0 

c initial  guess 

Z3a.«Sl 

Ot  L T A »FRRf)R 

IFCNC.r.T.n  Z3«Z-,01 
7 1 * Z 3-OEL  T A 
Z?«Z  3*r>tLTA 
VALl«FfZl) 
v ai.?«f  rz?) 

VAL3*F  f Z3) 

RR  OLLj*C73-7n/(Z2-Zn 
AL»H1 • f z *-7?1 / (Z?-Z1) 

&1 «ALAMI«*?»VAL 1-0FL1*a?«VAL?*(ALAM1*DEL1) *VALS 
CI*aLam1«(Al«m1*VAL1-.)El1*vaL2wVAlJ) 

NC?*'C2*1 

»•!  . 

TMGl.LT.O)  M-1. 

At A-i.(-?.,0EL1*vAl3) /(G1*5aS0RT  ( ABS  (61 •*?-«, *0ELl*Cl*VAL3)>) 


noon 
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V?«oOS  rtf**) 

Ct>  Tf 

30  Vt»TA\|H(EToH) 

ve*l  . 

3?  TfjP.Si*fT.wi*£?.*If2).V2 

( Z ) * F 2 * v 1 ♦ E T • V 2 

fiZ.5l»  t7»^)*'lTii-F.T»(El/E2)«COS(2*D)*(TOP) 

P E T o 8 N 

END 


FUNCTION  ETA(Z) 

COMMON  E?,F1,h,D,wavEL,PI,S1 
$ 1 « 1 . 

C IS  t^A*  IMAGINARY? 

IF  fH-r2-Z*»?.LE.tf1  Sl»-I. 
ETA .SCOT  f A«S(£I-E?-Z**a)I 
PETuPN 
End 


FUNCTION  XI(Z) 

COMMON  E?,Fi ,H(0,wAVEL,PI,SI 
T]=si-JKT  (El-l.-7**2) 
peTcpn 
eno 


Subroutine  sscz.et.axj 
Cj^.KJr.  E?,Ei,H,u,WAVEL,PI,Si 

Tnis  SjM(JiufjiJT  T NK  COMPUTES  K»".  SINCE  A«Y"«BE"  AND 
Fs*4rY»...HF»«  THIS  is  a tOUATION  IN  TWO 

UNK  IjO  S • r.JVfcN  F s 1 fxn  ANO  F $2  (X2)  FOR  SOME  VALUES 

of  (*y»»,e")  t he  slope  is  computed, 

C FS»FS(YY',ETA',Xl'»Kv"»E") 

Xl»FS(7,RT,AX,l.E-3,l.E-3) 

*2»fS(7,ET,AX,-1.E-3,1,E-3) 

Pm«  (Y  1 ♦*?)  /2.E-3 
»■ (XI -I ,E-3«RM)/1 .E-3 
slope \ ».pm/a 

C NOW  COMPUTE  THF  SLOPE  FOR  E (EQUIVALENT)  * * ’,VS.  l** 

SLOPE?*! .-ZaSlOPEI 

Type  sp.ei  ,z.slo°f:j  ,slopE2 

50  FuRmat  ( * » , Vt  1 « , tit  l . « , 2X  , flHWY '■  ,G1 1 ,« ,2*  , TMSLOPEl* , 

itil  !.4,2-<iTh$lOPF2»,GI  ».«) 
return 
EnO 


Function  fs(Bi  ,ii?,B3,H«,B5) 

C 31  ■KY*,B3»FTA',Hl8Xl',HaoKY",B5»E»' 

Common  R?,F 1 ,H,D,wAVEL,Pl,3l 
DI»3?P(h1 ,P2,6«,R5) 

D2»h3P  fM ,R3# a«f 85) 

C H?PiETA",H3P«XI»» 

Tip»3?*T»hh(H2*hT ♦ F 2 • B 3 

DjP»H3.E?»*?tTANH(H2*M)»E2*B2 

T2P»01*(TAMH(ti?»H)*b2AH«  ( (I  ,/COSM(B2«h)  )**2)  )pE2*0* 


39 


z»z3*»iAM.(^-z(n 

V*L'J«F(Z) 

CP!T.AHb( (7-7J)/Z3) 

IF  (Aps(VAl  .«)  .LE.FPS1  .0R.CRIT.1_E.EP3?)  GO  TO  98 

C CONVERGENCE  nt‘TAlNEO 


IF  (NC?.GT.IMAX)  GO  TO  97 

C IF  CONVERGENCE  FAILS  GO  to  97,  if  NEITHER  OF  the  ABOVE 

C rn»'i’ITTONS  ARE  SATISFIED  , obtain  NErt  POINTS  ANO  GO  back  (99) 

7l*ZP 
Z?*ZJ 

n»z 

VAl.l«V*L? 

VAl.  PbVAlI 
VaU3*VAL4 
GO  TO  99 


C CONVERGENCE  naS  FaILEO.TVPE  OUT  MESSAGE  ANO  TERMINATE  (1*1) 

97  TYPE  if» 

1*  FORMAT ( * *, 'ITERATION  OVER  LIMIT') 

TyRt  11 ,VAL«#C»IT 

11  F ! : R i'  / T ( * *,SHVAL4»*GU,<*,3X,5HCRIT»,G11,8) 

Type  1?. 7,73,7?, 71 

1?  FoRMAT('  '.?MZs,Gl«.7,?x,3HZ3«,Gl<*.7,PX,3HZ2»,Gl4,7,2X, 

UmZI  «,(;14.7) 

GO  TO  1/1 

c C 0 n v F F P E n C ( OtjTArrwEO.  COMPUTE  KT'* 

98  TYPt  1 J.VAI  u.CRIT 

13  f0r>aT('  ',5hvAL«»,G11.4,3K,5hCRIT»,G11,«) 

Tyre  ja.NCP 

1«  F0FTA1C  ',  1V1HITEKATI0N«,  14) 

ET«ETA  1 2 7 
AX*xI (Z) 

TYPf  is, 7, 8J ,ET, AX 

18  FORMAT  f ' ',3hky*,G14,7,3X,G9,2,2X,4hETA«,G14,7,2K» 

l 3mxI«,g1 4.7) 

C II.SFFT  PLOTTING  BLOCK  AND  OTHfR  OUTPUT  STATEMENT. 

C COMPOTF  EFFECTIVE  dielectric  constant 

EU*E 1 ”Z  a »2 
Typ»-  18.F1.E0 

18  Fup'ATf'  '.XHE1«,G12.5,3X,3HEU»,G12,3) 

c.ai  L SS(7,FT,  AX) 

c soi  vf  fof  slope  ok  kt" 


10J  CONTINUE 

1P1  STOP 

End 


Function  F(Z1 

rOMMON  C?.f 1 «H*0. HAVEL. PI. SI 
ET«ETA(Z) 

C IF  ETA'  Is  IMAGINARY  COSH  BECOMES  COS  ETC,, 

T F ( S I ) J1.3l.3P 
31  Vl«SlN(tT«H) 


I 


o?p«b i«f  »w»(  u ,/cnsH(«e*H) ) *a2)402*E2»»2 

1 »TAMM(t>?«R)*0l»F2 

f'>«<<i«fSTMfl,i*ni*,)i»n«roS(Rl*nn-t*?»(B'5»C0S(Pl»0) 

i-U»si.r.i.P)*c-R.,)*(Tip/niP)-ui*ti*cas(Bi*0)* 

1 CTlP/Ol»,)-H2»Kl»cni(Hl*D)*(fT?P/UJP)-{TlP*02P) 

1 / (f'lP««2)  I 
1?  fc  T U » N 
END 


FUNCTION  h?P  rR  1 , bP  , , R5) 

CjO,<o  j F->.F  1 ,M,n,WAVEL,PI,31 
oB<|)/(2,*B2) 

PE  T jpm 
E NO 


Function  n^pfai  ,bj,b*,P5) 
CtiRPOiv  E?,F1  .H,C,U4VEL,PI,  91 
B5P«(SS-2,«Ht  »«4)/(2,*B3) 
PtTUPN 
End 


